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Abstract — We propose an algorithm for identifying vertices 
from three dimensional (3D) meshes that are most important 
for a geometric shape creation. Extracting such a set of vertices 
from a 3D mesh is important in applications such as digital 
watermarking, but also as a component of optimization and 
triangulation. In the first step, the Ordered Statistics Vertex 
Extraction and Tracing Algorithm (OSVETA) estimates 
precisely the local curvature, and most important topological 
features of mesh geometry. Using the vertex geometric 
importance ranking, the algorithm traces and extracts a vector 
of vertices, ordered by decreasing index of importance. 

Index Terms — 3D shape, curvature estimation, discrete 
geometry, mesh optimization, vertex extraction. 

I. Introduction 

Curvature plays an important role in the global shape 
recognition, but there are others significant features of 
discrete object representations, which are involved in 
description of local topological characteristics. Thus, for 
example a low and high value of dihedral angles 
characterizes local topology ridge and furrow respectively. 
In the other hand, great distance of one point to its neighbors 
probably represents topological error, despite to evident 
high value of local curvature. Estimation of such features 
importance in combination with a local and global 
estimation of curvature is the main idea of this paper. In 
computers graphics 3D object is described in discrete form 
and fundamental element is a vertex given by its Euclidian 
coordinates. Therefore, the goal of our approach is a 
selection of important vertices, i.e. vertices that are most 
resistant to optimization and simplification processes. 

For a best defining criteria it is necessary a theoretical 
study of existing curvature estimation algorithms for the 
beginning. Results will demonstrate the effect of the 
different curvature estimation to the description of 
topological features. Within next phase follows an 
exploration of most important characteristics of optimization 
and simplification processes. In order to decrease a number 
of vertices, the optimization process uses topological criteria 
that determine "irrelevant" vertices for deletion from the 
mesh. We need to explain those useful criteria as principles 
of the vertex removal process, thus we define and 
successfully eliminate such vertices from our selection. 
Actually, in the context of vertex assessment well 
understanding of simplification criteria determines criteria 
that should be give bad grades. Thus, our approach provides 
forming criteria for action from both sides of problem. First 
set of criteria eliminates vertices from ordering based on 
matching their geometrical features with criteria conditions. 
The purpose of second set of criteria is ordering and finally. 



determination of the vertex extraction order in our approach 
is in direct proportion to its invariance in relation to the 
simplification and optimization process. 

The paper is organized as follows: Section II gives 
preliminaries with discussing a prior work in relevant areas 
of 3D geometry. In this section, we also define the notion of 
curvature and its sign, and introduce the concepts and 
terminology used in the paper. Typical methods of 
optimization and study of their characteristic properties are 
also located in this section. Section III presents the main 
idea as well as the details of OSVETA extraction algorithm. 
This section also provides conclusions about requirements 
that assessment criteria have to meet. In Section IV we show 
numerical results, which are obtained using our OSVETA 
software tools with real 3D meshes. In this section, we also 
discuss the limitations and possible shortcomings of the 
proposed approach and tradeoffs among complexity, 
distortion and security. Section V provides the summary of 
the numerical results, conclusions and future research 
directions. In Appendix (A and B), one may find statistic 
results of experimental test carried out to determine the 
rating of certain criterion. 

IT Preliminaries 

Estimation of curvature has long been unavoidable in the 
process of shapes perception and recognition [1], [2]. In 
1954, Attneave [1] published an important paper in which 
he observed that extrema of curvature along contours 
provided much of the information necessary to recognize 
objects from line drawings. Attneave manually picked 
points that corresponded to extrema of curvature and 
connected the points with straight lines. Remarkably, 
recognition based on such figures was found to be simple 
for human observers. However, defining important vertices, 
which form a simple and computer-recognizable 3D shape, 
is not trivial process. Thus, it is necessary to first define the 
problems and difficulties of discrete curvature estimation, 
but also point out differences between known methods. 

A. Discrete curvature estimation methods 

The notion of curvature originates from Gauss and his 
work in differential geometry [3] where he was first who 
formulated the main curvature features. Further curvature 
evaluation [4] methods improve mathematical interpretation 
and introduce new conditions, generalizing computation for 
all types of 3D meshes [5] [6]. For example, boundary 
meshes, that Gauss did not consider in his aforementioned 
work. However, the main discussions are conducted about 
accuracy of methods for discrete mesh computations. 
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1) Curvature Estimation using Differential Geometry 
From differential geometry [7] we know that for manifold 
surface M in ]R 3 , and each point on the given surface, one 
can locally approximate the surface by its tangent plane that 
is orthogonal to the normal vector n. For every unit 
direction e in the tangent plane, the normal curvature K is 
defined as the curvature of the curve that belongs to both the 
surface itself and the plane containing both n and e. The two 
principal curvatures: k x and k 2 of the surface M, with their 
associated orthogonal directions ei and e 2 are extreme values 
of all normal curvatures (Fig. 1. a). The mean curvature k h 
is defined as the average value of the normal curvatures 
/C//=(?C| +k 2 )/2. The Gaussian curvature k g is defined as the 
product of the two principle curvatures k g =k x k 2 . 




Figure 1. An illustration of (a) An infinitesimal neighborhood on a 
continuous surface patch; (b) 1-ring neighbors of point v,- and angles 
opposite to an edge; (c) External angles of a Voronoi region. 



For a given point v,- of discrete surface M the mean 
curvature normal (Laplace-Beltrami operator) 
K(v 1 -)=2/c / /(v 1 -)n(v 1 -) gives both the mean curvature 2 k h (x) 
and unit normal n(v,) at the vertex v,. Expressions for the 
mean curvature normal and the Gaussian curvature of a 
discrete surface that depend only on a vertex position and 
angles of adjacent triangles respectively are given in [8]: 
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where #f is a number of adjacent triangular faces at the point 
v„ and Oj is the angle of j - th adjacent triangle at the point v,-. 
A is the area of the first ring of triangles around the point 
v,. Minimizing a discretization error, Mayer, Desbrun, 
Schroder, and Barr [8] suggest a region A Mixed as a 
combination of Voronoi region for non-obtuse triangle and 
the half (or quarter) of a obtuse triangle region (Fig. 2.). 
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The expressions for discrete principal curvatures at the 
point v, are: 

V 0, ) = ( v,. ) + ^(v,.), k 2 ( v,. ) = k h ( v,. ) - TMv”) (4) 

where are 

A ( v i) = 4(v < )-* , g( v /). k h ( v /) = t|| K ( v /)||- (5) 





Figure 2. An example of mixed areas around a vertex use the midpoint of 
the edge opposite to the obtuse angle as the center: (a) 

Aon-Voronol = A and ( b > A on-Vorom,i =A> 2 - 



2) Fitting Quadric Curvature Estimation 

Fitting Quadric Curvature Estimation is one of the most 
used methods for curvature features computation today. 
Mclvor and Valkenburg give mathematical and theoretical 
background and comparison to previous similar methods in 
[9]. Later, Petitjean published survey [10] with comparison 
to the discrete differential geometry method, whereas 
Toolbox Graph (free Matlab software tools by Piere Alliez) 
is the practical application of this method; author refers to 
detailed algorithm in [11] and [12]. The idea of this method 
is that smooth surface geometry can be locally approximated 
with a quadratic polynomial surface. Method fits a quadric 
to point in a local neighborhood of each chosen point in a 
local coordinate frame. The frame is positioned at the 
chosen point and with the Z-coordinate axis, which is 
aligned along an estimated surface normal at that point. The 
curvature of the quadric at the chosen point is taken to be 
the estimation of curvature for the discrete surface. 

For a simple quadric form z' = ax' 2 +bx'y'+cy’ 2 procedure 
is given [9] by: 

• Estimation of surface normal n at the point v by one of 
two ways: simple or weighted averaging, or finding a 
least squares fitted plane to the point and its neighbors. 

■ Positioning of a local coordinate system ( x',y',z ' ) at the 
point v and aligning axis z' along the estimated normal. 
Mclvor and Valkenburg suggest aligning of the x' 
coordinate axis with a projection of the global x axis onto 
the tangent plane defined by n. This results in attitude 
matrix rotation 



y ~ [ r i,r 2 ,r 3 ] r (6) 

■ from the global coordinate frame to the local coordinate 
frame. Then: 
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where I is the identity matrix, and i is the global x axis 

[1,0, Of. 

■ Mapping the world data to the rotated principal frame 
with: x' = A'(x m , - v H- ) , where index w denotes that 

coordinates are expressed in global (world) coordinate 
frame. 

■ Fitting the mapped data to the rotated principal quadric 
z'= a'x' 2 +b'x'y'+c'y' 2 , and solve the resulting system of 
linear equations, giving a',b',c': 
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Computing both principal curvatures kj and k, using: 

k, = a' +c' +J(a' — c'Y +b ' 2 

, I , ,2 (9) 

k 2 = a! + c' - J(a' -c') +b'~ 

after that one can give expressions for Gaussian and mean 
curvature respectively: 

k g = Aa'c'-b' 2 , k h = a' + c' (10) 

Giving a new estimate of the surface normal at point v: 

<>■) 

if we use improvements suggested by Mclvor and 
Valkenburg and fit the mapped data to extended quadric: 
z = a'x 2 + b' xy + c' y 2 +d'x + e’y (12) 

Finally, we have estimation for Gaussian and mean 
curvature: 
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where a',b',c',d',e' are the parameters of the last quadric 
fitted. 

3) Comparison of curvature estimation methods 

A comparisons of mentioned methods that have 
performed by Mclvor and Valkenburg, and later Petitjean or 
Garimella and Swartz [13] gives enough information about 
accuracy. Our approach however does not require going into 
details. 

More purposeful is experimental computation of 
characteristic curvature features using both methods, and 
checking stability of vertices with highest calculated values 
of the all curvature features. For our considerations, the 
method is better than other one if its most important vertices 
survive the deletion process of optimization. 

Experimental and statistical results, presented in 
Appendix A, show different discrete values of characteristic 
curvature features for both methods (Table VI). Moreover, 
the difference of stable vertices appears, although only 
difference in their quantity is expected. 



B. Curvature estimation improvement 

Roughly speaking, the basic difference between 
mentioned methods is a direction of computation of 
characteristic features. Thus, fitting quadric method first 
calculates both principal curvatures, and then derives 
Gaussian and mean curvature. The differential geometry 
method however begins calculation with the Gaussian 
curvature. In discrete case this feature is calculated as a sum 
of all adjacent triangle angles at vertex, and then divided by 
an amount of the selected surface area around the chosen 
vertex. 

For our approach more important is the method with a 
direct impact on the calculation of Gaussian and mean 
curvature. Following our experimental results, we are going 
to improve the accuracy of curvature estimation for the 
method based on discrete differential geometry. We have 
assured that the local curvature computation using the 
Voronoi area did not give good experimental results; for 
obtuse angles, but also for any type of angles in strict local 
curvature estimation. Therefore, instead of Voronoi area we 
suggest using a barycentric area. 



Aarycenter X (cOt 0 + COt (X ) \\ X , - X f SUT 0 (14) 

6 yw,(/) 

Barycentric area at vertex v, is calculated as a surface 
inside a closed contour formed by the midpoints of all edges 
belonging to the vertex v, and the barycenters of its adjacent 
triangles. However, experimental results (See Appendix B) 
indicate a greater dependence of curvature estimation from 
the angle 6 than the theoretical region convergence. Indeed, 
this argument is valid for measuring the stability of certain 
areas only, not for the global curvature evaluation. 

C. Optimization and Simplification 

The second group of vertex features that we have used in 
consideration is a set of features that characterize salience of 
regions [14] [15] and resistance to transformations. Since we 
know that rotation, translation and uniform scaling of 3D 
models do not affect the topological features of mesh 
vertices, optimization and simplification most destructive 
affect the object geometry. Therefore, it is necessary to 
understand principles of optimization and simplification for 
a reliable vertex importance determination. The commonly 
used optimization tools include the algorithms such as: 
Mesh Optimization [16], Progressive Meshes [17], 
Simplification Envelopes [18], Quadric Error Metrics [19], 
Memoryless Simplification [20], and Quadric -based 
simplification [21]. These processes drastically change the 
geometric structure, removing vertices and deforming 
triangular faces. Comparisons of simplification algorithms 
are given in [22] and [23]. 

The basic principle of 3D surfaces simplification, that is 
also the most relevant to this paper, is the vertex decimation 
process. Schroeder, Zarge, and Lorensen [24] give the 
details of this algorithm in a prominent paper "Decimation 
of triangle meshes". The algorithm begins by multiple 
passes of testing the local 3D mesh geometry over all its 
vertices. During the test each vertex of 3D mesh is a 
candidate for removal, but vertex and its surrounding 
triangular faces are deleted if tested vertex meets specified 
decimation criteria. Resulting hole in the mesh is patched by 
a local triangulation. The vertex removal process is repeated 
until a termination condition is met, and steps of the 
algorithm are: 

■ Characterization of the local vertex, 

■ The decimation criteria evaluation 

■ The resulting hole triangulation. 

1) Local geometry characterization 

The first step of decimation algorithm is a 
characterization of local geometry and topology for a chosen 
vertex. The outcome of this process determines whether the 
vertex is a potential candidate for deletion. Each vertex can 
be classified according to one of three possible 
classifications: simple, complex and boundary vertex. The 
simple vertex can be further classified as a smooth vertex, 
an interior edge or a corner vertex. 

The simple vertex is surrounded by triangles where each 
edge containing the vertex is common to exactly two 
triangles. If the edge containing the vertex is shared with 
more than two triangles, or if vertex is common with at least 
one more than the cycle triangles, the vertex is complex. If 
an edge, which is branch of only one triangle, exists, then 
the vertex that belongs to this edge is the boundary vertex. 
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If the dihedral angle between two adjacent triangles is 
greater than a specified feature angle, then a feature edge 
exists. When a vertex is used by two feature edges, the 
vertex is an interior edge vertex. If three or more feature 
edges contain the same vertex, this common vertex is 
classified as the corner vertex. The illustration of 
characterization one may see in following figure [24] : 




(simple) (complex) (boundary (edge) (comer) 

Figure 3. Classification of vertices according to their local geometric and 
topological features. 

2) The decimation criteria evaluation 

The evaluation step determines which the vertices and 
triangular faces can be deleted, or eventually replaced by 
another triangulation. The main criterion for the decimation 
of a simple vertex is its distance from the average plane 
formed by adjacent vertices. If this distance is less than a 
given one, then the vertex is deleted (Fig. 4. a). Boundary 
vertex and the interior edge decimation use the distance 
from the edge criteria. In this case, the algorithm determines 
the distance from the line, defined by the two vertices that 
create the boundary or interior edge. If the distance is less 
than a certain value, vertex can be deleted (Fig. 4. b). 





Figure 4. Criteria of decimation: a) vertex distance from average plane, b) 
vertex distance from boundary 

Other simplification approaches are different from the 
specified vertex decimation approach in one, or even in all 
steps of algorithm. However, all of simplification algorithms 
minimize an energy function [25], and in this process some 
vertices, edges and faces are classified as risky primitives. 

3) Risky primitives and topological errors 
There is the list of risky geometric primitives: 

■ Isolated vertices, independent in space, 

■ vertices which belong to only one boundary edge, 

■ complex vertices and edges that join more than two faces, 

■ crossed edges without mutual vertex 

■ vertices with highest quantity of neighbors, 

■ vertices that are connected by collinear edges, 

■ boundary vertices, edges and faces, 

■ vertices, edges and faces in flat and smooth areas (d— >0). 
This listed group of primitives can be divided in two sub- 
groups: topological errors (the first four) and regular risky 
primitives (the last four). Topological errors are treated by 
our assessment criteria with strict zero grade, while regular 
risky primitives require different treatment and their 
assessment criteria should be adjustable to user. 

D. Curvature sign 

Based on the calculated Gaussian curvature, both 
principal curvatures and relationship between them, one can 
define vertex features in all areas of a given mesh M. The 
point v,e M is referred to as an umbilics if /ci(v,)=K 2 (v,'). If 



ki(x,)=k 2 (x,)=0, V,- is noted as a planar point. If k g (v,)=0, 
but v, is not a planar point, then v, is a parabolic point. If 
k g (v,-)> 0, v, is an elliptic point, and if k g (v,)<0, v, is called a 
hyperbolic point. Taking into account signs of both the 
Gaussian and mean curvatures we may distinguish the 
following eight basic surface shapes. 



TABLE I. Eight basic locally surface types in respect with their 
Gaussian and mean curvatures 
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Vertex features 


Surface shape 


3 £ 
li 'll 

O O 


planar 


planar 


o o 
k * 
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ridge convex X \ 


o o 
Si Si 
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valley concave 


o o 
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Si Si 
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Si Si 
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Si Si 


hyperbolic 


saddle 


o o 

Si Si 


hyperbolic 


saddle predominantly ridge 
convex 


o o 

v A 

o to 
Si Si 


hyperbolic 


saddle predominantly valley 
concave 



Shape of shown areas determines the condition of the 
vertex stability, whose measure is defined by a value of the 
curvature; more precisely, by characteristic features of 
curvature. 

III. Ordered Statistics Vertex Extraction and 
Tracing Algorithm (OSVETA) [28] 

The essence of OSVETA is represented by three steps: 
defining and ranking assessment criteria, accurate curvature 
evaluating and computing its characteristic features, as well 
as tracing the importance of extracted vertices in relation to 
the mesh topology. Based on the previous theoretical 
consideration, the algorithm first extracts all of the 
important geometric features. Statistical results of 
experimental tests are however crucial for criteria ranking, 
and finalizing a list of assessment criteria. Roughly 
speaking, the most important question is: which highest 
grades are the most significant? 

A. Assessment criteria ranking 

The following table gives the answer to the question and 
shows extracted criteria, and appropriate rates: 



TABLE II. OSVETA ASSESSMENT CRITERIA AND THEIR RATING 



N° 


Criterion 


Description 


Rate 


1 


Wmtn>= 0 


Positive minimal dihedral angle 


1.0 


2 


9< 2n 


Small Theta angle 


1.0 


3 


Kqi > 0 


Positive Gaussian curvature 


1.0 


4 


Vmax >= 0 


Positive maximal dihedral angle 


0.9 


5 


9>2n 


Big Theta angle 


0.8 


6 


k g < 0 


Negative Gaussian curvature 


0.8 


7 


Kgi < 0 


Negative Gaussian curvature 


0.7 


8 


k g > 0 


Positive Gaussian curvature 


0.4 
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yj max and y/ mi „ denote maximal and minimal dihedral 
angles at the vertex respectively. The angle 8 is the sum of 
all Ojj angles at the vertices (see Fig. 1. b). k g and k g , 
represent Gaussian curvatures estimated using respectively 
the both of method: discrete differential geometry method, 
and fitting quadrics method. The criteria ranking is 
performed using an automatic iterative computation of the 
survival vertices for the certain criteria selection. The results 
of tests are shown in Appendix A (See Table V). 

B. Irrelevant and Risky vertices elimination 

The elimination of some unimportant vertices from the 
assessment, and their sorting is also a significant process. 
Irrelevant and risky vertices are isolated from further 
computations and an outcome of elimination is decreasing 
of irrelevant vertices assessment probability. In addition, it 
reduces chance that topological errors to be considered in 
calculations; as it is mentioned previously, this kind of 
vertices (See Section II. C. 3) can get good grades for an 
eventual high value of curvature. 

Hence, in this step the algorithm extracts the matrix of 
topological error vertices E M , the matrix of boundary 
vertices B M , and a vector of adjacent vertices: N A for each 
vertex v,-. Using these matrices as well as the sign and a 
threshold value 1 of k g (x , ,) and k^x,) we first extract all risky 
vertices. The matrix of risky vertices is actually a union of 
matrices E M , B M 2 and a set of vertices which satisfy the 
condition: /c G (v,) = 0 a k/Xv,) = 0. 

In Table III we introduce new criteria with geometric 
features, obtained from conclusions of the consideration in 
Section II. C.: C e i is a coefficient of maximal elongated face; 
Cdug is a length of the longest connected edge to the vertex; 
Ctup is a value of the maximal angle of adjacent triangles; 
Cyis is a distance of the vertex from the average plane. V k g 
and V k h are gradients of Gaussian and mean curvature 
respectively. 



TABLE III. Irrelevant vertices elime'Jation criteria 



N° 


Criterion 


1 


(k g > 0 & k g < 7.05 13e-13) | (k g < 0 & k g > -2.8223e-9 ) 


2 


k h > 2.3688 | ( k h < 6.0769e-6 & k h > 0) 


3 


k h < -51929.7144 | ( K/f > -3.8486e-6 & k h < 0) 


4 


Kmax > > 39.5904 


5 


Kmax < -594.6627 1 (K raax < 0 & K max > -2.3077e- 15) 


6 


Xmin > 0.409451 (ttmin > 0 & K m in < 0.00017504) 


7 


x'min < -51843.002 | (K m i n < 0 & K m j n > -3.6061e-6) 


8 


8> 545.3552 | 0<213 | 0 = 360 


9 


kgi> 0.83844 \k gi < -1.2989 


10 


k H i > 1.2481 | ( Km >= 0 & k H i < 1 .671 e-5) 


11 


k hi < -2.3336 | (k H i < 0 & k h1 > -1.5326e-8) 


12 


Cel > 1.7682 I (Cel > 0 & Cel < 8.4903e-7 ) 


13 


Cel < -1.3404 1 (Cel < 0 & Cel > -6.9662e-5) 


14 


Cdug > 6.1607 | Cdug < 1.01 1 1 


15 


Ctup > 179.9596 | C TUP < 32.7365 


16 


Cv, s > 0.92554 | Cws < 3.7982e-6 


17 


V k g > 1467169.264 | V k g < 6.0536e-17 


18 


V k h > 2.8209| Vr„< 5.6712e-5 


19 


yfmax ' > ~ 90 | l//max < - 2 



1 Threshold values of all elimination criteria, including Gaussian and 
mean curvature values, are adjustable. In the OSVETA software [28], we 
allow three setting options: safe, extended and aggressive. 

2 Boundary vertices are often important for shape creation, and 
algorithm leaves to an user the choice of their elimination from ordering. 



C. The flowchart of Algorithm process 

The flowchart of all OSVETA algorithm processes is 
shown on Fig. 5.: 




M(VE) 




Figure 5. The flow chart of the OSVETA. 

For a mesh of a given 3D surface M(V,F) the result of 
OSVETA are two vectors: s, the vector of vertex stabilities 
arranged in a decreasing order, and i, the vector of 
corresponding indices. The mesh vertices are ordered with 
respect of decreasing stability form the vector V Q . V 0 =Vi, and 
the length-Z vector p is obtained by taking the first L 
elements from V 0 . N denotes the total number of vertices and 
also the total number of iterations for the shown process. 



IV. Numerical Results 



In this section we present a practical usability of our 
software in computation of the vertex importance, i.e. vertex 
invariance to the optimization process. We have used four 
3D mesh models (See Fig. 6.): 

(A) Christ the Redeemer [27], (B) Myron of Eleutherae [27], 
(C) Naissa by Bata [26], (D)Venus de Milo [27]. 

For brevity, we refer to these objects as to A, B, C, and D. 
All four models are different in total number of vertices and 
faces, but also in geometric structure. More precisely, they 
are different in percentage of curved and flatten areas. 
Myron of Eleutherae and Naissa by Bata are complex mesh 
models; both containing closed elements as subobjects. 
OSVETA algorithm works equally well with both 
homogeneous and complex meshes, as well as with open 
meshes with a boundary. 





Figure 6. Original 3D meshes: (A) Christ the Redeemer, (B) Myron of 
Eleutherae, (C) Naissa by Bata [18] and (D) Venus de Milo. 
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In order to compare OSVETA with a method of random 
vertex selection, we have used 'Optimize' modifier from 3D 
Studio Max 2012 application [29]. We have performed 
seven different levels of optimization with the following 
face thresholds: 

(A)(FT): 0,2,4,6,8,10,13,25 (B)(FT): 0,2,4,6,8,10,13,27 
(C)(FT): 0,2,4,6,8,10,13,26 (D)(FT): 0,2,4,6,8,10,13,24 

Higher maximal values are used for meshes with the 
larger total number of faces. The zero FT value does not 
affect the reduction of vertices and faces, and we consider it 
as an initial state of mesh geometry. The optimization with 
the maximum FT value, destroys completely the geometric 
structure of the mesh, leaving only 5-10% of the total 
number of vertices (See Fig. 7.). 




(A) (B) (C) (D) 

Figure 7. Destroyed geometries of 3D meshes after optimization with high 
level of FT. 



A determination of face threshold limits, preserving the 
usability of 3D mesh, is illustrated in Fig. 8. in the example 
of the Venus de Milo mesh. As can be seen from Fig. 8. the 
optimization face threshold higher than 8 leads to significant 
perceptual and geometric degradation of mesh. Thus we set 
the OSVETA criteria for assessing vertices without 
considering FT optimization values over 10. 




(00) (02) (04) 




(06) (08) (10) 

Figure 8. Perceptual degradation of the Venus de Milo 3D mesh as a 
function of the face threshold level. The subfigures correspond to the 
threshold levels 0,2,4,6,8,10. 

Experimental tests of stability with 1000 vertices selected 
by our algorithm showed the superiority of our approach 
compared to the randomly allocated group of 1000 vertices. 
The ratio of removed vertices from mentioned mesh models, 



Volume xx, Number x, 20xx 

denoted by (A), (B), (C) and (D) numerically is shown in the 
Table IV. Taking into account that different vertices are 
removed due to different levels of optimization, the table 
shows the accurate matching undeleted vertices of our 
algorithm. 



TABLE IV. Relation of number of deleted vertices by 

OPTIMIZATION, USING RANDOM AND OSVETA EXTRACTION ALGORITHMS. 



3D Model 


FT=0 


FT=2 


FT=4 


FT=6 


FT=8 


FT=10 


A# Vertices 


27802 


16627 


9789 


5991 


3743 


2432 


Random 


0 


502 


718 


838 


907 


948 


OSVETA 


0 


8 


51 


170 


333 


531 


B# Vertices 


100681 


49931 


32353 


21498 


15078 


11620 


Random 


0 


523 


686 


794 


848 


877 


OSVETA 


0 


6 


8 


19 


51 


94 


C# Vertices 


33465 


25334 


17472 


12146 


8588 


5870 


Random 


0 


276 


495 


654 


760 


845 


OSVETA 


0 


3 


21 


44 


77 


156 


D# Vertices 


17350 


12209 


6953 


3926 


2315 


1448 


Random 


0 


332 


622 


781 


872 


920 


OSVETA 


0 


1 


30 


147 


332 


522 



The one of most important features of our OSVETA 
extraction is stability of vertices grouped at the beginning of 
the selection. This fact is shown in the figure below used 
Naissa by Bata 3D mesh model, optimized using the limit 
level of face threshold (FT=10). Curve on the figure 
represents values of the Gaussian curvature at vertices 
selected using our algorithm (see Fig. 9.). 
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Figure 9. Gaussian curvature values at the first 50 extracted vertices by 
OSVETA: stabile vertices (filed circles) and deleted vertices by 
optimization (outlined circles). 

We observed another important characteristic from 
OSVETA computations: probability of two consecutive 
vertices deletion is very low. For example, for all 3D mesh 
models, using all optimization levels, in 1000 extracted 
vertices there is no consecutive deleted vertices. 

The presented vertex extraction approach based on 
OSVETA algorithm practically does not have restriction in 
use. A lack of such extraction is a dependence of curvature 
estimation accuracy, but we partially solve this problem 
using both evaluation methods as assessment criteria ( k g and 
k g\)- 

Although the computations of feature angles (0 and t//) are 
used in Gaussian and mean curvature estimation, 
experimental results of feature angles criteria indicate a 
fairly more efficiency (related on the extracted vertices 
stability) than Gaussian and mean curvature criteria (See 
Table V in Appendix A). 
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V. Conclusions and Future Work 

We have elaborated a novel method of extraction of 
important vertices from a 3D mesh. In several steps of the 
algorithm, the geometric importance of mesh vertices are 
ranked by the criteria based on the curvature and other noted 
important features. Thus, the two crucial goals are realized: 
first, the vector of irrelevant and risky vertices is separated 
from the mesh; second, the vector of the ordered vertex 
stability, with corresponding indices is extracted. 

Ordered Statistics Vertex Extraction and Tracing 
Algorithm have reached good efficiency showing invariance 
of extracted vertices to deletion by optimization. Flowever, 
we could see that a success of algorithm is determined by 
results of two research directions: curvature estimation and 
simplification of the process characteristics and effects. 
Thus, OSVETA future improvements will involve more 
precise mesh geometry estimation and better curvature and 
topological feature estimation. We believe that these 
enhancements would result in more accurate identification 
of stable vertices and consequently significant reduction of 
deletion probability. 

Appenddc A 

In this part of the paper, we present comparison of 
numerical results, which we have obtained from an 
experimental test. For the computation test, we have used 
four mesh models (See Section IV). 

The used models contain a different total number of 
vertices and faces. Geometric structure of each objects is 
also unique, i.e. the ratio of curved and flat areas is distinct 
for each model. In addition, Myron of Eleutherae and Naissa 
by Bata are complex mesh models, i.e. the both contain 
closed elements as subobjects. 



TABLE V. Criteria efficiency for given number of selected 

VERTICES. 



N° 


Criterion 


Number of selected vertices 


1 


2 


3 


4 


5 


10 


50 


too 


500 


1000 


1 


K C > 0 


1 


2 


2 


3 


4 


7 


19 


27 


86 


148 


2 


k g < 0 


0 


1 


2 


3 


4 


7 


21 


31 


86 


128 


3 


Kh> o 


1 


2 


2 


2 


3 


5 


13 


23 


37 


37 


4 


Kh< o 


1 


1 


1 


1 


1 


2 


4 


9 


34 


69 


5 


Kmax > 0 


1 


2 


2 


2 


2 


4 


17 


25 


86 


129 


6 


^-min ^ 0 


0 


0 


0 


1 


1 


1 


4 


9 


42 


80 


7 


6 > 360 


0 


1 


1 


2 


2 


7 


31 


54 


138 


192 


8 


0 <360 


1 


1 


2 


3 


4 


7 


31 


47 


130 


185 


9 


k G i > 0 


1 


2 


3 


3 


3 


5 


21 


38 


99 


131 


10 


Kqi < 0 


0 


0 


0 


0 


1 


3 


17 


31 


82 


125 


11 


k hi > 0 


0 


1 


1 


1 


1 


2 


11 


18 


79 


116 


12 


Km < 0 


0 


1 


1 


2 


2 


2 


12 


20 


68 


99 


13 


Cel > 0 


0 


0 


0 


0 


1 


1 


6 


10 


41 


70 


14 


Cel <= 0 


0 


1 


1 


1 


1 


1 


2 


3 


24 


34 


15 


Cvis > 0 


0 


0 


0 


1 


1 


4 


8 


18 


52 


87 


16 


Ctup > 0 


0 


0 


0 


0 


0 


1 


4 


9 


39 


71 


17 


Cdug > 0 


0 


1 


1 


1 


1 


1 


1 


3 


18 


34 


18 


V Kr > mean 


0 


1 


1 


1 


1 


1 


9 


14 


45 


86 


19 


V Kh < mean 


0 


0 


0 


0 


0 


0 


0 


0 


7 


19 


20 


V Kq > mean 


1 


2 


2 


3 


3 


5 


13 


27 


84 


151 


21 


V Kg < mean 


0 


0 


0 


0 


0 


0 


0 


0 


0 


4 


22 


¥4n a x>=0 


0 


0 


1 


1 


1 


2 


11 


18 


112 


204 


23 


V^min ' > 0 


1 


2 


3 


4 


5 


6 


22 


37 


86 


134 



First, we test an efficiency of all criteria, i.e. we select and 
sort vertices of source mesh by 23 criteria and find for each 
criterion a number of selected vertices, matched to vertices 
of optimized object. The test is repeated for a different 
number of selected vertices. Previous table (Table V) shows 
results of the test for Venus de Milo 3D mesh object 
(FT=13). 

In the next experimental test, we compare selections of all 
noted criteria. Each of 23 criteria selects ten vertices and 
Fig. 10. represents the histogram of frequent selected 
vertices. 




Figure 10. Frequency of the vertex selection (Venus de Milo). 

As we can see, only two same vertices are selected by 3 
criteria, three same vertices are selected by two criteria, and 
rests of vertices are selected by unique criterion. 

The Histogram also indicates that the most of criteria 
selects unique vertices. Moreover, some similar criteria, e.g. 
Gaussian curvatures, estimated by two mentioned methods, 
select the complete different first 10 vertices. Table VI 
shows vertices selected by k g > 0 and k gi > 0 criteria for the 
maximal level of optimization (FT = 13). 



TABLE VI. Differences of selected vertices by similar criteria. 



N° 


Criterion 


Selected vertices 


1 


K C > 0 


8755 


8769 


8773 


9386 


9439 


10374 


10376 


9 


Kgi > 0 


5827 


5829 


6174 


6833 


15925 







Indeed, geometric position of selected vertices in 3D 
space is the most important information. Thus, 3D plot on 
Fig. 11. shows both sets of selected vertices in their different 
topological regions. Perceptual and topological importance 
of regions to whom belong selected vertices is also clearly 
shown of next figure. We note that there is no identically 
selected vertex by both of criterion, i.e. vertices selected by 
each criterion are at different topological regions. 




(a) (b) (c) 



Figure 11. 3D position of first 50 vertices selected by Kg and Kqi criteria, 
(a) Kq > 0 - circles, (b) zoom of (a), (c), Kgi > 0 - diamonds. 
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Appendix B 

We compute for Section II. B extraction efficiency of k g > 
0 and k g < 0 criteria, using Voronoi and Barycentric regions 
in curvature estimation. For a better comparison with 
previous results Venus de Milo is used in extraction. 
Numbers of extracted vertices using both compared regions 
are given in following table: 



TABLE VII. Criteria efficiency for given number of selected 

VERTICES. 



N° 


Criterion 


Number of selected vertices 


1 2 3 4 5 10 50 100 500 1000 


Curvature estimation using barycentric area in calculation 


i 


k g > 0 


1 


2 


2 


3 


4 


7 


19 


27 


86 


148 


2 


Kq< 0 


0 


1 


2 


3 


4 


7 


21 


31 


86 


128 


Curvature estimation using Voronoi area in calculation 


i 


Kq> 0 


1 


2 


2 


3 


4 


6 


18 


25 


84 


140 


2 


Kq< 0 


1 


2 


3 


3 


4 


4 


12 


20 


80 


123 



Significant differences in criteria efficiency using 
barycentric and Voronoi area in curvature estimations are 
observable only with large numbers of selected vertices. 
This fact is sufficient argument for use barycentric area in 
Gaussian and mean curvature calculation, because use of a 
small number of selected vertices are very rare. 

In addition we may conclude that the use of barycentric 
areas in curvature calculations gives better results in the 
hyperbolic i.e. saddle shapes of mesh surfaces (k g < 0). In 
elliptic and parabolic regions (k g > 0 and k g = 0 
respectively) curvature calculations, using both types of 
areas, ensure results with similar accuracy. 
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